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Precisely engineered tunnel junctions exhibit a long sought effect that occurs when the energy 
of the electron is comparable to the potential energy of the tunneling barrier. The resistance of 
metal-insulator-metal tunnel junctions oscillates with an applied voltage when electrons that tunnel 
directly into the barrier's conduction band interfere upon reflection at the classical turning points: 
the insulator-metal interface, and the dynamic point where the incident electron energy equals the 
potential barrier inside the insulator. A model of tunneling between free electron bands using the 
exact solution of the Schrodinger equation for a trapezoidal tunnel barrier qualitatively agrees with 
experiment. 



PACS numbers: 85.75.-d, 72.25.-b, 75.70.-i 

Tunneling through a barrier is one of the most fun- 
damental problems in physics with profound technolog- 
ical implications [H HJ [3J 0] |S]. Electron tunneling can 
be realized in multilayer systems consisting of two con- 
ducting electrodes separated by an insulating material. 
These tunneling devices can have electrodes that are nor- 
mal metals {e.g., Au), superconductors (Al), or ferro- 
magnets (Fe), while the insulators range from semicon- 
ductors (Ge) to metal oxides (MgO). In all cases, the 
energy difference between the Fermi energy of the elec- 
trodes and the conduction band of the barrier defines 
the height of the tunneling barrier. The relative en- 
ergy of the tunneling electrons is varied by applying a 
voltage bias between the electrodes. The barrier mate- 
rial's properties generally play no role because the bar- 
rier potential is typically much greater than the energy 
of the tunneling electron. However, recent technologi- 
cal advances have led to robust oxide barriers that can 
withstand large electric fields, which allows access to elec- 
trons with energies comparable to the barrier potential. 
It was predicted long ago that a tunneling electron could 
access electronic states in the barrier and undergo res- 
onant tunneling [(F, though this has not been directly 
observed in metal-insulator- metal junctions with a single 
tunnel barrier. This Letter reports well defined and re- 
producible bias-dependent oscillations of the differential 
resistance in CoFeB/MgO/NiFe tunnel junctions that are 
consistent with interfering electrons within the MgO bar- 
rier. Further, we are able to use the ferromagnetism of 
the electrodes to investigate these oscillations in terms of 
a tunneling magnetoresistance. Qualitative agreement 
with these data is obtained using a simple tunneling 
model where the electrodes are treated as free electron 
bands and the tunneling matrix elements are determined 
by solving the Schrodinger equation exactly for a trape- 



zoidal barrier. 

The procedure used to fabricate our magnetic tunnel 
junctions (MTJs) was reported previously [7J. We in- 
vestigated about thirty MTJs each of CoFeB/MgO/NiFe 
(NiFe) and two sets of CoFeB/MgO/CoFeB (CFB1 and 
CFB2), for a total of nearly ninety junctions. As de- 
posited CoFeB is amorphous, while NiFe is polycrys- 
talline. The barriers were formed by oxidizing 16 A of 
Mg, which should produce a 13 A thick MgO barrier. All 
devices were fabricated identically with the exceptions of 
the free layer material and the procedure for oxidation 
and annealing. While data presented here are for 1 /im 2 
MTJs, our observations were independent of device shape 
and area, which varied from 1100 nmx420 nm ellipses to 
600 nm— 10 /im diameter circles; heating and series resis- 
tances were thus negligible. The temperature dependence 
can be attributed to thermal smearing [8], and all devices 
satisfied the MT J tunneling criteria [9l [10] , proving that 
tunneling is the primary conduction mechanism. 

Direct measurements of the differential resistance 
(dV/dl) were made with the magnetizations of the fer- 
romagnets (in remanence) in parallel (p) and antiparallel 
(ap) configurations using a high resolution ac resistance 
bridge and standard lock-in techniques. The dc bias was 
applied to the free layer (NiFe or CoFeB) with the pinned 
layer (CoFeB) grounded. Fig. [T] shows dV/dl measure- 
ments of a NiFe device at 5 K. An obvious oscillation for 
positive biases is apparent in the ap state; a low ampli- 
tude oscillation is also present in the p state. The dV/dl 
evolve continuously between these states as a function of 
the angle between the free and pinned magnetizations. 
dV/dlp equals dV/dI ap at several biases, the highest of 
which (+1.8 V) was observable at 5K because of an in- 
crease in the breakdown voltage at low temperatures. 
These crossing biases were independent of temperature 
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FIG. 1: (Color online) Differential resistance measurements 
of NiFe devices at 5 K in the p and ap magnetic configura- 
tions, (inset) The bias dependence of dI p /dV for NiFe (green 
squares) is fit well by a 4 th order polynomial (black line). 



and the relative orientation of the ferromagnets. 

Though the underlying physics is contained in the 
dV/dl for each magnetic orientation, inspecting the dif- 
ferential junction magnetoresistance (dMR) is a conve- 
nient way to emphasize this oscillatory behavior. Defin- 
ing the dMR as 



dMR 



dV/dI ap - dV/dI p 
dV/dl 



it is most common to observe dMR > because the den- 
sity of states bottleneck typically causes dV/dI ap to ex- 
ceed dV/dI p . Thus, the most striking feature of Fig. [2| 
is that the dMR of the NiFe devices oscillates about 
zero when electrons tunnel from CoFeB into NiFe. In 
contrast, no oscillations were observed for either CFB1 
or CFB2 (though these did have dMR < at high bi- 
ases). The symmetry of the biases for which dMR ini- 
tially changes sign (i.e., the first zero-crossings for both 
polarities) mimics that of the barrier shape. The barrier 
parameters extracted using the Brinkman, Dynes, and 
Rowell (BDR) model [TT] for each device type are shown 
schematically in Fig. [2] CFB1 had a relatively symmet- 
ric barrier and symmetric zero-crossings. The first zero- 
crossing for NiFe occurs when electrons tunnel toward 
the low barrier interface. CFB2 had the opposite sym- 
metry of the NiFe devices both in barrier heights and 
zero-crossings. Based on this seemingly general behav- 
ior, we predict similar oscillations would be seen for both 
polarities in NiFe/MgO/NiFe devices. 

One possible explanation for this behavior is dynamic 
resonant tunneling mediated by electron interference, a 
result of Fowler-Nordheim (FN) tunneling [T5]. In the 
FN regime, the incident electron energy exceeds the po- 
tential near the collector within the barrier, which im- 
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FIG. 2: (Color online) Normalized differential magnetore- 
sistance for NiFe (bottom, blue), CFB1 (middle, red), and 
CFB2 (top, black) at 300 K. The nominal dMR(V=0) and 
resistance-area products were 80% and 3.1 kf2 • /im 2 for NiFe, 
and 120% and 6.0 kQ-fim 2 for both CFB devices. The barrier 
heights (eV) from BDR fits of dI p /dV are indicated schemat- 
ically to the right, and a cartoon of electron standing waves 
is shown to the left. 



plies that electrons tunnel directly into the MgO con- 
duction band (see cartoon in Fig. [2]). These electrons 
can thus be treated as plane waves near the collector 
interface. Incident and reflected electrons may then in- 
terfere and establish standing waves within this region 
of the barrier (where the kinetic energy of the electron 
is real). The maximum oscillation amplitude exists in 
the ap state, which is the result of spin-dependent re- 
flection from the MgO-NiFe interface. Spin-up electrons 
tunneling from the CoFeB emitter with positive magne- 
tization toward the NiFe collector with negative mag- 
netization are preferentially reflected at the MgO-NiFe 
interface because of a spin bottleneck in the density of 
states. When the applied bias allows for resonant tunnel- 
ing in the ap state, spin-up emitter electrons dominate 
the tunneling current. The result of this is that the an- 
tiparallel state conductance exceeds that of the parallel 
state, which leads directly to dMR < 0. Negative dMR 
indicates that tunneling is dominated by minority spin 
electrons rather than majority electrons (as defined in the 
collector electrode) , and oscillations suggest that the rel- 
ative conductance of the spin species changes with bias. 
Gundlach showed that oscillations of the conductance (in 
tunneling between normal metals) can only be described 
by exactly solving the Schrodinger equation, and that 
this is evidence for the failure of the simplest application 
of the Wcntzcl-Kramcrs-Brillouin (WKB) approximation 
[6]. This effect has been identified in semiconductor sys- 
tems |13j and scanning tunneling microscopy |14j . but 
not in metal-insulator-metal heterostructures. 

Figure [3] shows that these experimental results are 
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FIG. 3: (Color online) Model results for the dMR bias depen- 
dence for NiFe (blue line), CFB1 (red dots), and CFB2 (black 
dashes) show qualitative agreement with experiment. For a 
thickness of 28 A and no roughness, the barrier heights (eV) 
used in the model are indicated schematically to the right. 



qualitatively reproduced by a model that uses the ex- 
act solution of the Schrodinger equation with spin-split 
free electron bands representing the ferromagnetic elec- 
trodes. A free electron model is justified here because 
tunneling is dominated by s-like electrons [Tj)| EH IT?] . 
This model reproduces the angular dependence of the 
experimental dMR when the polarization of either elec- 
trode is varied as cos 9, including a nodal behavior of 
the biases where dMR equals zero (i.e., when dV/dI p = 
dV/dI Qp ). Asymmetric barriers were approximated as 
single-thickness trapezoids with barrier heights <p\ and 
03 at the pinned and free layer interfaces. We included 
material-dependent effective masses, and thus imple- 
mented the BcnDanicl-Dukc boundary conditions: con- 
tinuity of the wave function ipi(xij) — ipj(xij), and 
flux d x ipi(xi.j)/m* — d x ipj(xi,j) /m* , at Xij between 
materials i and j |18j . Because ab initio calculations are 
unavailable for disordered electrodes, we approximate the 
(spin- independent) effective masses (m* = 1.3rn e ) and 
band bottoms (2.3 eV below the Fermi energy, Ep) for 
both electrodes by Fe band structure calculations [T5] , 
The MgO effective mass was taken to be 0.4m e [20]. The 
spin splittings (A) of the electrodes were estimated via 
Ai = E F [(l - Pi) 2 /{I + Pi) 2 - 1], where i denotes the 
material, and the polarizations (P) of CoFeB and NiFe 
were taken to be 55% and 45% from tunneling spin po- 
larization measurements |211 122] . The tunneling currents 
in the parallel and antiparallcl states were calculated in- 
cluding the spin-dependent densities of states |23] via 

l p cxj (p^N^Nf +D u n{ 1} n1 3) ) F(E)dE, 

I ap (x J (D u iY T (1) 7Vf ) +D lt N^N^)F(E)dE, 



where subscripts denote the pertinent spin sub-bands, su- 
perscripts denote the materials, and F(E) is the differ- 
ence in Fermi functions of the two electrodes [f^(E) — 
f( 3 \E-eV)}. The spin-dependent densities of states 
for materials 1 (grounded, pinned layer) and 3 (bi- 
ased free layer) are respectively N$ — Nm\E) and 
Nn — Nn (E — eV), where subscripts represent the 
spin sub-bands. The tunneling matrix elements D mn — 
D mn (s, 0i , 03, V, E) are the probabilities of tunneling be- 
tween spin sub-bands m and n in electrodes 1 and 3, 
respectively, and were obtained by exactly solving the 
Schrodinger equation. Spin-flip processes were neglected. 
The upper limit of the transverse integration was trun- 
cated at 2% of the Fermi energy because the tunneling 
current is dominated by wave vectors within a cone of 
~ 8° from normal incidence 24J. This model is simple 
and therefore appealing, though more intricate analyses 
|25] may capture more effectively the underlying spin- 
dependent physics. 

The location and amplitude of the negative dMR peak 
shift only slightly when interfacial roughness is included 
in the calculation in a previously demonstrated manner 
[2"rj] . This is because the bias of the initial oscillation 
is set by the barrier height at the collector interface (the 
threshold bias for FN tunneling equals this height) . Note 
that for CFB1 or CFB2 oscillations were not observed 
because the barrier heights were comparable to the max- 
imum applied bias. On the other hand, a more significant 
effect is seen for the positive peak around 1.4 V. Includ- 
ing 15% roughness, which was the roughness determined 
for similar MTJs [27], causes the amplitude of the latter 
peak to fall from 58% to 45%, and requires a mean thick- 
ness of 35 A to keep the dMR peaks at biases consistent 
with the data of Fig. [2] the discrepancy from the growth 
thickness is reasonable for this qualitative model. 

The model needed to explain the present data implies 
barrier parameters that are different from those obtained 
by BDR fits. To qualitatively reproduce the data, the 
model requires similar 0i with only 03 significantly dif- 
ferent between the three MTJ types (shown schemati- 
cally in Fig. [3]). This is reasonable since the fabrication 
process for all the devices was identical until the barrier 
oxidation; perturbations due to this step should predom- 
inately affect 03, not <\>\. The BDR fits, on the other 
hand, indicate a roughly constant average barrier height 
with both 4>i and 03 different for each MTJ type, which 
is less reasonable considering the fabrication procedure. 
The model calculations use a thickness of 28 A (when 
neglecting roughness), while the BDR fits yield ~ 8 A. 
This discrepancy is most likely due to interfacial rough- 
ness and the WKB approximation failing in these devices 
(as indicated by the electron interference presented here). 

FN tunneling is required for the interference of elec- 
trons within the barrier. In this regime, both the width 
and average height of the barrier at the Fermi energy de- 
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crease with increasing bias. The result of this should be 
an exponentially increasing conductance for applied bi- 
ases greater than the collector interface barrier height. 
The parallel state conductance data are fit well by a 
fourth-order polynomial (Fig.[I]inset), showing that devi- 
ations from the low-bias parabolic behavior exist, but the 
data are not exponentially increasing at the highest bi- 
ases. The delayed transition to exponential conductance 
may be related to the effective mass of the tunneling elec- 
tron, interfacial roughness, or may be an emerging char- 
acteristic of coherent tunneling through crystalline bar- 
riers, possibly originating from heating effects [271 155] , 

Similar oscillations are obtained in resonant tunneling 
studies where a normal metal spacer exists between the 
barrier and one ferromagnet (see e.g., Ref. [21J] and refer- 
ences therein). We refer to those as "static" phenomena 
because the thickness of the normal metal is fixed. The 
phenomenon we report here is "dynamic" because the 
thickness of the interference region can be tuned by the 
applied bias: electrons inside the MgO conduction band 
can be treated as free electrons, making this region di- 
rectly analogous to the normal metal spacer used in static 
studies. Additionally, the present case allows the effec- 
tive mass of the MgO conduction band to be estimated 
from the oscillation period. Recalling that the deBroglie 
wavelength (in A) is A = 0.529em e /m* , where e is the 
dielectric constant (e ~ 3) and m*/m e is the reduced 
mass of MgO [30], the measured oscillation period cor- 
responds to A ~ 2.6 A, and m* ~ 0.6m e , in reasonable 
agreement with the expected 0.4m e [20 . 

Alternate explanations for our observations can be 
ruled out. The oscillations are odd functions of bias, and 
thus cannot be explained by emission phenomena such as 
magnons |31j . The persistence of the oscillatory nature 
at finite bias with interface roughness excludes localized 
barrier states and interface resonant states |32[ 133]. The 
angular dependence is strong evidence that this effect is 
not due to quantum size effects in the electrodes [54"] . 
A similar negative dMR region in Co2MnSi/MgO/CoFe 
tunnel junctions was recently interpreted as an energy 
gap in the Co2MnSi minority-spin band [35] . but this 
explanation does not apply to our case because no such 
gap is expected for NiFe. While it is not possible to rule 
out unknown density of states effects, dynamic resonant 
tunneling is the most convincing origin because of its sim- 
plicity and ability to reproduce numerous experimental 
features from different MTJs with different materials. 

In summary, novel bias-dependent oscillations were ob- 
served in the differential resistances of mctal-insulator- 
mctal (CoFcB/MgO/NiFe) tunnel junctions. These long 
sought oscillations are due to dynamic resonant tunneling 
mediated by interference of electrons that tunnel into the 
conduction band of the insulator (MgO). A coherent tun- 
neling model using the exact solution of the Schrodinger 
equation and free electrons representing the electrodes 
qualitatively reproduced the bias dependence. The tun- 



ability of this newly demonstrated phenomenon, as well 
as its spin-dependence, may help advance the develop- 
ment of tunable resonant tunneling systems for funda- 
mental spintronics physics and applications. 
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